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ABSTRACT 

The  atomic  structures  of  several  Ca-Mg-Cu  and  Ca-Mg-Zn  amorphous  alloys  have  been 
detennined  using  ab  initio  molecular  dynamics  simulation  and  neutron  diffraction.  Partial  pair 
distribution  functions  have  been  produced  and  the  pair  bond  distances  and  partial  coordination 
numbers  have  been  reported  for  these  alloys.  Similarities  and  differences  in  the  amorphous 
structures  of  the  Ca-Mg-Cu  and  Ca-Mg-Zn  alloys  have  been  discussed. 
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1.  INTRODUCTION 

Ca-Mg-based  bulk  metallic  glasses  (BMGs)  have  unique  properties.  They  are  built  on  two  simple 
metals,  Ca  and  Mg,  while  all  other  BMGs  are  transition-metal  based  alloys  [1,2].  Many  Ca-Mg- 
based  alloys  have  very  good  glass  fonning  ability  (GFA),  and  even  their  ternary  alloys,  e.g.  Ca- 
Mg-Zn  or  Ca-Mg-Cu,  can  be  made  fully  amorphous  when  they  are  cast  in  ~8-10  mm  thick  plates 
[3, 4, 5, 6, 7, 8]  or  up  to  15  mm  diameter  rods  [9, 10].  These  glasses  have  very  low  Young’s  and 
shear  moduli  that  are  comparable  with  the  moduli  of  human  bones  [11,12],  low  density  ( 1 .6-2.8 
g/cm  )  [1],  and  strong  relaxation  dynamics  of  the  super-cooled  liquid  [13].  Glass  fonning  ability, 
as  well  as  physical  and  mechanical  properties  of  metallic  glasses,  is  believed  to  depend  on  the 
type  of  short  range  order  (SRO)  and  medium  range  order  (MRO)  of  alloying  elements  in  the 
amorphous  structure.  It  is  therefore  interesting  to  analyze  the  amorphous  structure  of  Ca-Mg- 
based  BMGs. 

Structural  analysis  of  metal-metal  BMGs  is  mainly  focused  on  transition  metals  glasses  such  as 
Zr-Cu  [  14, 15, 16, 17],  Zr-Pt  [18],  and  Zr-Cu-Al  [19,20].  The  results  indicate  that  the  packing  of 
atoms  in  these  materials  is  not  random,  but  is  strongly  influenced  by  chemical  interactions.  In 
particular,  icosahedral  SRO  has  been  identified  and  correlated  to  their  good  GFA  [21].  On  the 
other  hand,  the  atomic  structure  of  Ca-Mg-Zn  BMGs,  identified  with  the  use  of  X-ray  and 
neutron  diffraction  and  Reverse  Monte  Carlo  (RMC)  simulation,  shows  no  icosahedral  SRO 
[22].  Instead,  five-fold  bonds  in  the  fonn  of  pentagonal  bi-pyramids  have  been  found  to  be  the 
most  populous  structural  units  in  these  BMGs  [22].  A  large  fraction  of  five-fold  bonds  and  the 
lack  of  icosahedral  SRO  has  also  been  found  in  the  molecular  dynamic  (MD)  simulated 
amorphous  structure  of  Mg-Cu  alloys  [23]. 

In  the  present  work,  the  atomic  structures  of  several  Ca-Mg-Cu  and  Ca-Mg-Zn  BMGs  were 
simulated  with  the  use  of  ab  initio  (quantum)  molecular  dynamics  (QMD)  simulation.  The 
simulated  structures  were  then  used  to  calculate  partial  (PRDF)  and  total  (RDF)  radial 
distribution  functions,  pair  bond  distances  and  partial  and  total  coordination  numbers.  The  QMD- 
simulated  results  were  validated  through  the  experimentally  detennined  RDFs. 

2.  EXPERIMENTAL  PROCEDURES 

QMD  simulation  of  the  amorphous  structures  of  the  Ca6oMgi5Zn25,  Ca6oMg25Zni5, 

Ca6oMgi5Cu25,  Ca6oMg25Cui5,  Ca5oMg25Cu25,  and  Ca40Mg25Cu35  alloys  was  conducted  using  the 
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Vienna  ab  initio  Simulation  Package  (VASP)  [24].  The  Projector  Augmented- Wave  (PAW) 
method  [25,26]  and  Perdew-Burke-Emzerhof  (PBE)  functional  [27]  were  used  to  describe 
interacting  valence  electrons.  Cubic  simulation  boxes  with  periodic  boundary  conditions 
contained  200  atoms,  and  the  box  size  was  detennined  by  the  density  of  the  modeled  material. 
The  simulation  was  conducted  at  T  point  only,  with  a  first  order  Methfessel-Paxton  smearing 
function  (c=0.2  eV).  The  liquid  was  equilibrated  at  1000K,  then  was  quenched  to  300K  in  100K 
temperature  steps,  equilibrating  at  each  temperature  for  3000  steps.  Another  3000  consecutive 
configurations  were  further  sampled  at  300K  with  a  time  step  of  5  fs,  and  averaged  PRDFs, 
gijQMD(r),  were  calculated  from  these  configurations  using  equation  (1),  thus  taking  into  account 
thennal  vibrations  and  assuming  no  structural  relaxation  within  15  ps. 


gy(r)  = 


1  dnyjr) 
\Tw1cjp0  dr 


(1) 


Here  d«ij  are  the  number  of  elements  of  type  j  between  distances  r  and  r+dr  from  an  element  of 
type  i,  cj  is  the  atomic  fraction  of  the  element  of  type  j,  and  pQ  is  the  alloy  density.  The  QMD 
models  provided  a  complete  set  of  six  simulated  PRDFs  for  each  of  the  studied  alloys,  which 
well  agreed  with  the  experimental  diffraction  data.  Due  to  the  limited  box  size  and  the  periodic 
boundary  conditions,  gijQMD(r)  were  calculated  at  r  <  9  A. 

Fully  amorphous  samples  of  the  studied  Ca-Mg-Cu  and  Ca-Mg-Zn  alloys  were  prepared  by  melt¬ 
spinning  in  the  fonn  of  ribbons  as  described  in  [22].  The  densities,  pQ,  of  the  amorphous  alloys 
were  measured  with  a  helium  pycnometer  AccuPyc  1330  VI. 03  and  the  values  (in  g/cnr  and 
atoms/A3)  are  given  in  Table  1.  Neutron  diffraction  experiments  were  conducted  at  room 
temperature  under  vacuum  using  the  General  Materials  (GEM)  diffractometer  at  the  ISIS  high- 
intensity  pulsed  neutron  source  (Rutherford  Appleton  Laboratory,  Didcot,  UK)  [28]  and  yielded 
the  neutron  total-scattering  structure  factors  (SSF),  S(Q): 

1  da  3 


S(Q)  = 


N  dO. 


“ZfA2 


(2) 


1=1 


Here - is  the  differential  neutron  cross-section  per  unit  solid  angle  Q  for  the  3-  component 

N  dQ 

alloy,  Cj  and  b,  are,  respectively,  the  atomic  fraction  and  the  coherent  bound  neutron  scattering 
length  of  element  i.  Each  S(Q )  was  Fourier  transformed  to  give  a  real  space  neutron  total  RDF, 
G(r),  defined  as  [29,30]: 
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n 


The  Lorch  modification  function  [31],  L(Q),  with  a  maximum  momentum  transfer  Qmax  of  25  A"1 
was  used  to  reduce  termination  ripples  in  the  Fourier  transform  (Equation  2). 

Using  Equation  (1),  the  partial  coordination  number  of  an  element  j  in  the  first  coordination  shell 
of  the  element  i  was  determined: 

finax 

Ny  =  4 ncjpo  |  g/- dr  (5) 

r  ■ 

'min 

Here  rmin  and  rmax  are  the  positions  of  the  start  and  end  of  the  first  peak  in  the  respective  gi/r). 
The  simulated  structures  were  statistically  analyzed  using  pair  and  three-body  correlation 
functions,  Voronoi  tessellation  and  nearest  neighbor  approaches,  which  allowed  calculation  of 
local  structural  features  (SRO  and  MRO)  such  as  total  and  partial  CN,  type  and  distribution  of 
characteristic  coordination  polyhedra,  atomic  volume,  and  local  packing  fraction. 

3.  RESULTS  AND  DISCUSSION 

3.1.  Structure  Factors  and  Radial  Distribution  Functions 

The  experimental  SSFs,  S(Q),  for  the  six  amorphous  alloys  are  given  in  Figure  1.  The  S(Q) 
curves  show  a  pronounced  peak  in  the  Q  range  from  2  to  3  A’1  and  several  small  diffuse  peaks  at 
higher  Q  values.  The  intensities  of  these  peaks  decrease  quickly  with  Q  and  no  oscillations  are 
apparent  beyond  ~15  A"1.  As  the  Cu  and  Zn  concentrations  increase,  the  first  diffraction  peak 
broadens  (Figure  1).  The  presence  of  a  pre-peak  at  Q  ~  0.8  -1.5  A’1  may  indicate  the  presence  of 
MRO  in  these  alloys  [32,33],  as  well  as  the  fluctuation  in  the  atomic  scattering  cross-sections: 
The  solutes,  Cu,  Zn  and  Mg,  have  larger  neutron  scattering  cross-sections  than  the  solvent,  Ca. 
The  experimental  RDFs,  G(r),  are  given  in  Figure  2  as  solid  lines.  Systematic  changes  of  the 
shape  of  the  first  RDF  peak,  which  is  located  at  r  ~  2. 1  -  4.5  A,  with  the  alloy  composition  are 
observed  (Figure  2).  For  example,  the  first  RDF  peak  in  Ca6oMgi5Cu25  has  a  maximum  at  r  = 

3.05  A  and  two  shoulders,  one  at  r  =  2.58  A  and  another  at  r  =  3.70  A.  As  the  amount  of  Ca 
decreases  and  the  amounts  of  Cu  and  Mg  increase,  the  intensity  of  the  first  shoulder  increases 
and  the  intensity  of  the  second  shoulder  decreases.  On  the  other  hand,  an  increase  in  the 
concentration  of  Mg  at  the  expense  of  Cu  results  in  the  intensity  of  the  second  shoulder  to  exceed 
the  intensity  of  the  former  maximum  (compare  RDFs  of  Ca6oMgi5Cu25  and  Ca6oMg25Cui5  in 


G(r)  =  X 

i,j= 1 


cicibibj[g.(r)-l]  = 


2n rpt 


J  QS  (Q)L(Q)  s\n(rQ)dQ 


(3) 
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Figure  2).  Replacing  Cu  with  Zn  noticeably  decreases  the  intensity  of  the  first  shoulder.  Because 
the  first  RDF  peak  is  a  sum  of  the  concentration-  and  scattering  length-  weighted  peaks  of  the  six 
relevant  PRDFs  (Equation  3),  such  a  strong  concentration  dependence  of  the  shape  of  the  first 
RDF  peak  is  apparently  due  to  different  partial  coordination  numbers  and  pair  bond  distances  in 
the  studied  alloys.  The  very  different  neutron  scattering  length  of  Cu  (/y  u  =  7.718  fin)  and  Zn 
(. bzn  =  5.680  fim)  need  also  to  be  taken  into  account  when  comparing  Ca-Mg-Cu  and  Ca-Mg-Zn. 
Figure  3  shows  six  QMD-simulated  PRDFs  for  the  CaeoMgisZms,  and  Ca6oMgi5Cu25  amorphous 
alloys.  The  PRDFs  are  rather  smooth  since  they  account  for  atom  vibrations.  The  first  peak  of 
each  of  the  PRDFs  corresponds  to  the  atom  pair  interactions  within  the  first  coordination  shell 
and  it  provides  information  on  the  atom  pair  bond  distance  distributions  and  partial  coordination 
numbers  around  the  alloying  elements.  Comparison  of  the  PRDFs  of  the  CaeoMgisZms,  and 
Ca60Mgi5Cu25  amorphous  alloys  shows  that  the  intensities  of  the  Ca-Ca,  Ca-Mg  and  Ca-Zn  peaks 
in  the  Ca6oMgi5Zn25  alloy  are  similar  to  the  intensities  of  respective  Ca-Ca,  Ca-Mg  and  Ca-Cu 
peaks  in  the  Ca6oMgi5Cu25  alloy.  At  the  same  time  the  intensities  of  Zn-Zn  and  Mg-Zn  peaks  are 
smaller  and  the  intensity  of  the  Mg-Mg  peak  is  larger  in  Caf^Mgi  sZrbs  than  the  respective 
intensities  of  Cu-Cu,  Mg-Cu  and  Mg-Mg  peaks  in  Ca6oMgi5Cu25-  This  interesting  observation 
may  indicate  that  the  Ca-Ca,  Ca-Mg  and  Ca-Zn(Cu)  interactions  are  similar  in  these  alloys,  while 
the  tendencies  for  the  Cu-Cu  and  Mg-Cu  bonding  are  stronger  than  for  the  Zn-Zn  and  Mg-Zn 
bonding.  Mg  atoms  seem  to  have  higher  tendency  to  cluster  in  the  Ca6oMgisZn25  alloy  than  in  the 
Ca60Mgi5Cu25  alloy. 

Figure  4  illustrates  the  composition  dependence  of  two  PRDFs,  gcu-Cu(r)  and  gca-Mg(r)  for  the 
four  Ca-Mg-Cu  amorphous  alloys.  The  effect  of  the  composition  on  the  position  of  the  first 
PRDF  peak  is  weak  relative  to  the  effect  on  the  positions  of  the  higher  order  peaks.  This 
indicates,  to  the  first  order,  that  bond  distances  between  respective  atom  pairs  in  the  first 
coordination  shell  are  not  affected  by  the  composition  in  the  Ca-Mg-Cu  system.  A  stronger 
composition  dependence  of  the  intensity  of  the  first  peak  for  the  Cu-Cu  PRDF  than  for  the  Ca- 
Mg  PRDF  may  indicate  that  the  Cu-Cu  chemical  order  is  stronger  than  the  Ca-Mg  chemical 
order. 

Using  the  QMD-simulated  PRDFs  and  Equation  (3),  simulated  total  RDFs  were  obtained  for 
every  alloy  and  compared  with  the  respective  experimental  neutron  RDFs  (Figure  2).  A  very 
good  match  is  seen.  Since  the  QMD  simulations  were  conducted  completely  independently  from 
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the  neutron  diffraction  experiments,  these  good  fits  indicate  that  the  simulated  PRDFs  are 
statistically  consistent  with  the  amorphous  structures  of  the  studied  alloys  and  they  can  be  used 
to  obtain  statistical  infonnation  on  the  effect  of  the  alloy  composition  on  the  SRO  in  these  alloy 
systems. 

3.2  Atom  Pair  Bond  Distances 

Using  the  QMD-simulated  PRDFs,  the  most  frequent  (mode,  r"ode ),  weighted  average  (mean, 
f.Mean  j  anc)  cut-off  (maximum,  r^u,~OJ}  )  distances  between  atom  pairs  in  the  first  coordination 

shell  were  calculated  for  the  studied  alloys  and  their  values  are  given  in  Table  2,  Table  3,  and 
Table  4,  respectively.  The  weighted  average  distances,  r™ean ,  between  the  i  and  j  atoms  were 

calculated  using  the  following  equation  [22]: 

rCut-off  I  Cut-off 

rMea"  _  J rjgij{r)dr  \gv{r)dr  (6) 

o  /  o 

The  rMode  values  (Table  2)  are  always  smaller  than  the  respective  rMean  values  (Table  3).  This 
indicates  that  the  first-shell  gij(r)  peaks  are  asymmetric  (non-Gaussian).  The  difference  is  larger 
for  the  Zn-  than  for  the  Cu-  containing  pairs,  which  may  indicate  a  stronger  interaction  of  Cu, 
than  Zn,  with  other  elements.  The  rMode  and  rMean  values  can  be  compared  with  atomic 
separations  in  crystalline  metals  and  alloys.  Crystalline  metallic  separations,  tm,  and  covalent 
distances,  rc,  were  estimated  as  a  sum  of  metallic  and  covalent  radii  of  respective  elements 
(Table  5).  The  metallic  radii  are  half  the  nearest  distance  between  atoms  in  crystal  lattices  of  the 
respective  pure  metals  [43].  The  covalent  atomic  separations  were  obtained  from  [44],  based  on 
the  analysis  of  the  nearest  distances  in  crystalline  intermetallic  compounds  reported  in  the 
Cambridge  Structural  Database.  From  Table  5,  the  rMean  and  rMode  values  for  Ca-Cu,  Ca-Zn,  and 
Mg-Cu  atom  pairs  are  smaller  than  ym-  However,  for  all  other  atomic  pairs  rMean  values  are  larger 
and  r Mode  values  are  smaller  than  rM-  At  the  same  time,  rMode  for  the  atom  pairs  containing  Cu  are 
even  smaller  than  the  respective  covalent  bond  distances  in  crystalline  compounds.  The  Mg-Cu 
(2.73  A),  Ca-Cu  (3.03  A)  and  Ca-Zn  (3.14  A)  mode  pair  distances  and  are  much  shorter  than 
those  calculated  from  the  Ca-Ca,  Mg-Mg,  Cu-Cu  and  Zn-Zn  rMode  values  (i.e.  2.82  A,  3.13  A, 
and  3.20  A,  respectively,  see  Table  5).  Shortening  of  Mg-Cu  and  Y-Cu  bonds  has  also  been 
observed  in  amorphous  Mg60Cu3oYi0  [34].  Such  bond  shortening  often  indicates  strong  chemical 
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interactions  between  these  atom  pairs.  No  bond  shortening  is  observed  between  Ca-Mg  and  Mg- 
Zn  pairs. 

The  observation  that  some  nearest-neighbor  bond  lengths  are  shorter  than  those  in  the  competing 
long-range  ordered  crystals  may  help  explain  why  these  Ca-Mg-Cu  and  Ca-Mg-Zn  alloys  are 
easy  glass  formers.  The  equilibrium  crystalline  phases  that  compete  with  the  amorphous  phase 
appear  to  minimize  total  energy  by  sacrificing  optimal  SRO  to  achieve  beneficial  long-range 
order.  On  the  other  hand,  the  absence  of  long-range  atomic  order  allows  metallic  glasses  to 
reduce  the  energy  difference  between  the  metastable  glass  and  equilibrium  crystal(s)  by 
optimizing  short-range  atomic  interactions  and  arranging  the  atoms  in  efficiently  packed  clusters 
[35,36,  37,38],  The  largest  contribution  to  condensed  phase  stability  is  expected  to  come  from 
nearest-neighbor  interactions  and  volume  minimization,  so  that  the  optimized  short-range 
interactions,  as  well  as  the  presence  of  MRO,  in  metallic  glasses  can  produce  structures  that  have 
only  a  small  energetic  disadvantage  relative  to  the  equilibrium  crystalline  structure  [21,22]. 
Kinetic  constraints  from  quenching  and  packing  frustration  restrict  the  long  range  atomic 
redistribution  needed  to  achieve  long-range  order  and  to  further  minimize  the  total  system 
energy,  thus  favoring  glass  formation.  The  present  results  suggest  that  metallic  glasses  not  only 
have  SRO,  but  may  have  “better”  combination  of  short-range  topological  and  chemical  order  (in 
terms  of  optimal  bond  length  and/or  atomic  arrangement)  than  the  competing  crystals. 

3.3.  Coordination  Numbers  and  Chemical  Short  Range  Order 

The  average  partial  and  total  coordination  numbers  around  Ca,  Mg  and  TM  atoms  (TM  is  Cu  or 
Zn),  from  the  QMD  simulations,  are  given  in  Table  6.  There  is  almost  no  effect  from  replacing 
Cu  with  Zn  on  the  partial  and  total  coordination  numbers  around  Ca  and  Mg  atoms  in  the 
Ca6oMgi5TM25  and  CaeoMg^sTMis  amorphous  alloys  (see  the  first  4  columns  in  Table  6).  The 
total  coordination  number  around  Ca  and  Mg  are  CNca  =  14.3  and  CNMg  =  1 1.3-1 1.5  for 
Ca60Mg25TMi5  and  CNCa  =  14.4  and  CNMg  =  1 1.9-12.0  for  Ca6oMgi5TM25.  The  number  of  TM 
and  Ca  increase  and  Mg  decrease  in  the  first  coordination  shell  of  Ca  and  Mg  clusters  with  an 
increase  in  the  concentration  of  the  TM  element  in  these  alloys.  At  the  same  time,  Cu  has  a 
slightly  different  neighbor  environment  than  Zn.  In  particular,  the  amount  of  Ca  is  slightly  lower 
and  Mg  slightly  higher  around  Cu  than  around  Zn.  This  is  in  agreement  with  the  earlier  discussed 
observation  that  the  Mg-Cu  interactions  are  stronger  than  Mg-Zn  interactions. 
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Increasing  the  Cu  concentration  from  15  to  35  at%  continuously  increases  the  total  coordination 
numbers  CNca  from  14.3  to  16.2  and  CNMg  from  1 1.3  to  12.8,  due  to  a  faster  increase  in  the 
number  of  smaller  Cu  atoms  (CNca-cu  and  CNMg-cu  increase  from  1.6  to  4.9  and  from  1.1  to  3.2, 
respectively)  and  a  slower  decrease  in  the  amount  of  the  larger  Ca  atoms  (CNca-ca  decreases  from 
9.2  to  7.0  and  CNMg-ca  decreases  from  8.5  to  6.9)  (Table  6).  On  the  other  hand,  CNCu  almost  does 
not  depend  on  the  alloy  composition.  However,  similar  to  the  environment  around  Ca  and  Mg, 
CNcu-Cu  increases  from  0.6  to  1.6  and  CNcu-ca  decreases  from  6.4  to  5.6  with  the  increase  in  Cu 
concentration.  The  number  of  Mg  atoms  around  Ca,  Mg  and  Cu  has  a  tendency  to  weakly 
increase  with  an  increase  in  the  concentration  of  Cu  (see  Table  6). 

To  detennine  if  the  atoms  in  the  amorphous  structure  are  distributed  randomly  or  they  have  a 
tendency  to  form  specific  clusters,  the  degree  of  chemical  short  range  order  (CSRO)  in  the  first 
coordination  shell  of  Ca,  Mg  and  TM  atoms  was  estimated  using  a  CSRO  parameter,  ai(jk),  which 
was  defined  for  ternary  alloys  in  ref.  [39]: 

am  =  1  -  (CNjj+  CNik)/ ((cj+Ck)CNi)  for  i  +  j  +  k  (8) 

Here  cj  is  the  atomic  fraction  of  the  element  j,  CNjj  is  the  partial  coordination  number  of  element 
j  in  the  first  coordination  shell  of  the  element  /,  and  CA)  is  the  total  coordination  number  of 
element  /.  Negative  a^k)  values  indicate  the  presence  of  CSRO,  i.e.  the  total  concentration  of  j 
and  k  atoms  in  the  first  coordination  shell  exceed  their  total  average  concentration  in  the  alloy. 
The  positive  values  indicate  the  excess  of  /-elements  in  the  first  coordination  shell  of  i  atom,  i.e. 
chemical  short  range  clustering  (CSRC)  of  atoms  /. 

The  a,i(jk)  values  for  the  studied  Ca-Mg-TM  glasses  are  given  in  Table  7.  The  parameter  ai(jk)  is 
positive  for  Ca  atoms  and  is  negative  for  Mg  and  TM  atoms.  This  result  indicates  that  Ca  atoms 
tend  to  cluster  with  each  other,  while  Mg  and  TM  atoms  have  a  tendency  to  attract  unlike  atoms 
in  the  first  coordination  shell  and  fonn  CSRO.  This  tendency  for  CSRO  is  stronger  for  the  TM 
elements  than  for  Mg  in  all  studied  alloys,  except  Ca6oMg25Cuis.  The  CSRO  of  the  TM  element 
with  Ca  and  Mg  increases  with  an  increase  in  the  concentration  of  the  TM  element  in  the  alloy. 
The  presence  of  CSRO  near  TM  and  Mg  and  CSRC  near  Ca  atoms  indicates  that  the  amorphous 
structure  of  the  Ca-Mg-TM  alloys  is  not  a  random  solution  of  the  alloying  elements,  but  it 
consists  of  TM-  and  Mg-  centered  atomic  clusters,  with  specific  CSROs.  Stronger  CSRO  around 
the  TM  elements  may  indicate  that  the  TM-centered  clusters  dominate  over  the  Mg-centered 
clusters. 
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3.4  Voronoi  Analysis  of  the  Atomic  Clusters  in  Ca-Mg-TM  Amorphous  Alloys 

Voronoi  tessellation  and  the  resulting  Voronoi  cells  allow  identification  of  characteristic 
coordination  polyhedra  present  in  an  amorphous  structure  [22,37].  A  coordination  polyhedron  is 
defined  as  an  /'-centered  cluster  with  vertices  at  the  first-shell  atom  positions  and  edges 
coinciding  with  the  interatomic  bonds  in  the  first  shell  [40],  Each  coordination  polyhedron  can 
be  assigned  a  Voronoi  signature  ( ns,n4,n,5,ri6 ),  where  nm  is  the  number  of  vertices  common  to  m 
polyhedron  faces  (or  edges)  [41].  m  is  also  called  the  vertex  coordination.  Those  /-centered 
clusters  with  the  same  Voronoi  signature  are  considered  to  be  topologically  equivalent  (even 
though  they  may  not  be  identical). 

Table  8  lists  Ca,  Mg  and  TM-centered  clusters  identified  in  the  amorphous  structures  of  the 
studied  Ca-Mg-TM  alloys.  Only  clusters,  which  are  present  at  a  fraction,^,  of  0. 1  or  higher  in  at 
least  one  of  the  alloys  are  tabulated.  Just  three  types  of  the  Ca-centered  polyhedra  have  the 
fractions  fP>  0.1.  These  are  (0,2, 8, 4),  (0,1,10,4)  and  (0,2, 8, 5).  The  total  fraction  of  these  clusters 
varies  from  0.36  in  Ca6oMg25Zni5  to  0.075  in  Ca40Mg25Cu35  and  many  other  different  types  of 
Ca-centered  polyhedra  are  present  at  /7  ~  0.05  or  less.  Many  different  types  of  clusters,  with  low 
fraction  of  each  of  the  type,  may  indicate  that  the  Ca-centered  clusters  are  not  likely  building 
blocks  for  the  amorphous  structures. 

Six  characteristic  Mg-centered  and  five  TM-centered  clusters  with  fp>0.\  are  also  shown  in 
Table  8.  All  of  them  except  (1,3, 3, 3)  are  Kasper-type  polyhedra,  which  account  for  36  to  80% 
among  the  Mg-centered  clusters  and  44  to  82%  among  the  Cu-centered  clusters  (these  clusters 
are  the  least  populated  in  Ca40Mg25Cu35  and  most  populated  in  Ca6oMg55Zn25).  The  high  fraction 
of  the  selected  TM-  and  Mg-centered  Kasper-type  polyhedra  should  indicate  that  the  amorphous 
structure  of  the  Ca-Mg-TM  alloys  mainly  consists  of  these  clusters,  which  is  in  agreement  with 
the  recently  proposed  ECP  model  of  the  amorphous  structure  [38,42].  Since  Kasper  polyhedra 
are  polytetrahedral,  and  many  other  polyhedra  also  contain  tetrahedra,  we  conclude  that  the 
topological  SRO  of  Ca-Mg-TM  is  polytetrahedral  in  nature,  which  has  been  demonstrated  in 
many  other  metallic  glasses  [22],  The  current  study  therefore  further  supports  the  idea  that  the 
SRO  of  MGs  is  characterized  by  polytetrahedral  packing  via  Kasper  polyhedra  and  their 
distorted  variants. 


9 


Approved  for  public  release;  distribution  unlimited 


The  analysis  of  the  most  populated  clusters  given  in  Table  8  indicates  that  five-coordinated 
vertices  dominate  in  all  the  clusters;  however,  the  amorphous  structures  of  the  Ca-Mg-TM 
glasses  contain  a  low  fraction  of  icosahedral  SRO  (Kasper  polyhedra  with  (0,0,12,0)  Voronoi 
index).  This  is  presumably  due  to  the  large  mismatch  of  atomic  sizes  in  the  nearest  neighbor 
shell  (e.g.,  rCa=1.97  A  and  rCu=1.28  A),  which  renders  the  distorted  Kasper  polyhedra  more 
topologically  favored.  A  large  fraction  of  live-fold  bonds  and  the  lack  of  icosahedral  SRO  have 
also  been  found  in  the  MD-simulated  amorphous  structure  of  Mg-Cu  alloys  [23].  The  dominance 
of  live-fold  vertices  suggests  that  tetrahedra  prefer  to  cluster  into  pentagonal  bi-pyramids. 

Table  8  also  demonstrates  the  different  local  topology  in  Ca6oMgi5Zn25  and  Ca6oMgi5Cu25 
(columns  1  and  3),  which  differs  by  only  the  type  of  the  TM  element.  Specifically,  because  Cu  is 
smaller  in  size,  and  has  obvious  bond  shortening  with  both  Ca  and  Mg,  it  is  expected  to  induce 
higher  CN  around  Ca/Mg,  as  well  as  more  distortion  in  the  Kasper  polyhedra.  Indeed,  the  most 
popular  cluster  type  around  Ca  is  (0,1,10,4)  with  CN=15  in  Ca6oMgi5Cu25,  and  (0,2, 8, 4)  with 
CN=14  in  Ca6oMgi5Zn25.  For  Mg  centered  clusters,  it  is  (0,3, 6, 4)  with  CN=13  in  Ca60Mgi5Cu25, 
and  (0,2,8, 1)  with  CN=1 1  and  (0,0,12,0)  with  CN=12  in  Ca6oMgi5Zn25. 

3.5  Local  volume  and  packing  fraction 

Assuming  atoms  are  hard  spheres,  the  global  packing  fraction  can  be  easily  calculated  using  the 
atomic  radii  and  the  density  of  the  alloy.  This  approach  yields  a  packing  fraction  of  0.76-0.80 
for  the  alloys  in  this  study  (see  Table  9).  However,  this  treatment  is  very  crude,  and  it  does  not 
provide  information  on  local  packing  fraction  specific  to  an  atom,  an  element,  or  a  cluster.  In  the 
following  we  propose  a  method  that  can  be  used  to  describe  the  local  volume  and  packing 
fraction. 

The  weighted  Voronoi  tessellation  assigns  each  atom  with  a  Vornoi  cell,  and  its  volume  can  be 
considered  as  the  atomic  volume  (i.e.,  the  volume  occupied  by  the  center  atom  at  its  position). 
The  summation  of  the  atomic  volume  defined  as  such  equals  to  the  total  volume  of  the  alloy. 
Therefore,  the  ratio  between  the  hard  sphere  volume  (calculated  from  the  metallic  radius)  and  the 
atomic  volume  (Voronoi  cell  volume)  is  on  average  smaller  than  1,  and  reflects  the  packing 
density  around  the  specific  atom  (i.e.,  atomic  packing  fraction  in  Table  9).  The  challenge  now  is 
how  to  define  the  packing  fraction  for  a  short-range  cluster,  and  the  key  is  to  find  out  a  way  to 
quantify  the  actual  volume  occupied  by  a  cluster  in  the  amorphous  matrix.  An  /-centered  cluster, 
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3 

consists  of  the  atom  i  in  the  center  and  CN .  =  V  CN,,  atoms  in  the  first  shell.  The  atoms  in  the 

i  ij 

j= i 

first  shell  are  shared  with  other  neighbor  clusters.  The  average  volume  of  the  cluster,  Vu  is  then 
defined  as: 


V,  =  vt 


/  +St 

7=  1  1 


2,  OV„V, 


+  K 


(9) 


Here  v,-  and  vy  are  the  average  Voronoi  volumes  of  atoms  /  and  j,  and  Ky  is  the  average  number  of 
clusters,  which  are  the  neighbors  to  the  /-centered  cluster  and  which  share  the  atom  j  in  their  first 
shell.  The  averaging  is  conducted  over  the  whole  number  of  clusters  of  type  /  presented  in  the 
amorphous  structure  and  vertices  of  type  j  in  the  cluster  /.  It  is  expected  that  Ky  =  1  +  jUy,  where 
Hy  is  the  average  coordination  of  vertex  j  in  the  /-centered  cluster: 


m= 3 

and  thus 


(10) 


v  v+f^lL 

7*2  +  11* 


(ID 


Here  fm  is  the  fraction  of  j  vertices  with  the  vertex  coordination  my  in  all  /-centered  clusters. 


Indeed,  in  general  two  neighbor  clusters  can  share  a  face,  an  edge,  or  a  vertex.  The  average 
number  of  faces  or  edges,  which  have  a  common  vertex  j  is  equal,  by  definition,  to  /Uy.  A 
geometric  analysis  shows  that  the  neighbor  clusters  having  a  common  vertex  j,  can  fill  space 
only  if  they  share y-vertex-containing  faces  with  the  /-centered  cluster  and  thus  the  number  of 
these  neighbors  should  be  equal  to  juy.  One  additional  neighbor  cluster  will  share  only  the  vertex 
j  with  the  /-centered  cluster.  The  juy  values  can  be  directly  obtained  from  the  results  of  the 
Voronoi  tessellation  of  the  amorphous  structure. 

Equations  (10)  and  (11)  were  used  to  calculate  the  average  vertex  coordinations  and  average 
volumes  of  the  Ca-,  Mg-  and  TM-  centered  clusters  for  all  the  studied  alloys.  The  results  are 
given  in  Table  10.  The  average  coordination  of  the  Ca-,  Mg-,  and  TM-  vertices  is  ~  5.1,  5.0,  and 
4.8,  respectively,  and  it  weakly  depends  on  the  alloy  composition.  At  the  same  time,  the  average 
cluster  volumes  are  noticeably  composition-dependent.  As  it  is  expected  from  the  numbers  and 
radii  of  the  atoms  in  the  clusters,  the  Ca-centered  clusters  have  the  largest  volume  (—97- 1 09  A3) 
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and  the  TM-centered  clusters  have  the  smallest  volume  (—53-64  A3).  The  Mg-centered  clusters 
have  an  intennediate  volume  of -74-84  A3.  The  cluster  volumes  in  the  Zn-containing  alloys  are 
larger  than  in  the  respective  Cu-containing  alloys,  probably  because  the  atomic  radius  of  Zn  is 
larger  than  that  of  Cu,  and  Cu  incurs  more  significant  bond  shortening.  In  the  Ca-Mg-Cu  alloys, 
the  volume  of  the  clusters  decreases  with  an  increase  in  the  concentration  of  Cu  and/or  decrease 
in  the  concentration  of  Ca/Mg  due  to  an  increasing  number  of  Cu  and  decreasing  number  of 
Ca/Mg  in  the  first  shell  of  all  three  type  clusters. 

The  average  packing  fraction  of  the  /-type  clusters  was  estimated  by  calculating  the  total  volume 
of  the  atoms  belonging  to  the  averaged  /-cluster  (using  the  equation  similar  to  (1 1),  where  the 
Voronoi  volumes  are  replaced  with  respective  atom  volumes)  and  dividing  it  by  V°ver .  Atomic 
radii  reported  in  Table  5  were  used  to  calculate  the  atom  volumes.  The  results  are  given  in  Table 
10.  The  Ca-centered  clusters  are  packed  slightly  more  densely  than  the  Mg-  and  Cu-  centered 
clusters.  In  Ca-Mg-Zn  alloys,  the  average  packing  fraction  of  Ca-centered  clusters  is  -0.77  and  it 
decreases  to  -0.74-0.75  for  Mg-centered  and  to  -0.74-0.76  for  Zn-centered  clusters.  The  packing 
fraction  of  the  clusters  is  slightly  smaller  in  the  Zn-containing  alloys  than  in  the  respective  Cu- 
containing  alloys,  probably  due  to  stronger  interactions  of  Cu  with  other  elements  leading  to 
more  intense  bond  shortening.  An  increase  in  the  concentration  of  Cu  in  the  Ca-Mg-Cu  alloys 
increases  the  clusters’  packing  fraction,  thus  supporting  the  presence  of  strong  Cu-Cu,  Cu-Mg 
and  Cu-Ca  bonding  in  these  amorphous  alloys. 

3.6  Electronic  structure  of  Ca-Mg-TM  amorphous  alloys 

It  has  been  shown  in  the  above  discussion  that  the  presence  of  Cu  or  Zn  in  the  Ca-Mg-TM 
amorphous  alloys  leads  to  different  bonding  characters  and  thus  properties.  It  is  therefore 
intriguing  to  examine  the  electronic  structure  of  Ca-Mg-Cu  and  Ca-Mg-Zn,  especially  the 
differences  associated  with  the  TM  element.  To  this  end,  site-  and  orbital-projected  electronic 
density  of  states  (DOS)  was  calculated  using  the  configurations  obtained  in  QMD  simulation,  on 
a  3x3^3  Monkhorst-Pack  grid.  Only  valence  electrons  (2  for  Ca,  2  for  Mg,  1 1  for  Cu,  and  12  for 
Zn)  are  accounted  for  as  the  low-lying  core  electrons  are  unlikely  to  participate  in  the  bonding. 
Wigner  radii  for  the  decomposition  are  1.947  A  for  Ca,  1.524  A  for  Mg,  1.270  A  for  Zn,  and 
1.312  A  for  Cu.  The  results  are  compared  in  Figure  5. 
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The  most  prominent  difference  between  Figure  5a  (Ca6oMgi5Zn25)  and  Figure  5b  (Ca6oMgi5Cu25) 
is  the  very  different  position  of  the  narrow  and  intense  3d-band:  for  Zn  it  is  over  7  eV  below  the 
Fermi  level,  while  for  Cu  it  is  2-3  eV  below  the  Fenni  level.  As  a  result,  the  3d-band  of  Zn 
barely  overlaps  with  the  s-  and  p-  conduction/bonding  bands  in  Ca/Mg,  whereas  for  Cu  the 
overlap  is  significant.  Such  overlap  between  the  3d-band  of  TM  and  the  sp-band  of  simple  metals 
is  known  to  result  in  strong  chemical  interactions  and  orbital  hybridization,  as  previously 
demonstrated  in  Al-Cu,  where  the  Al-Cu  bonds  are  shortened  by  over  6%  [19].  The  DOS  in 
Figure  5  may  explain  why  Cu  induces  more  obvious  bond  shortening  than  Zn  (e.g.,  Ca-Cu  bond 
is  shortened  by  7.26%,  while  Ca-Zn  bond  is  shortened  by  only  5.41%). 

CONCLUSIONS 

1 .  The  atomic  structures  of  six  Ca-Mg-TM  ternary  metallic  glasses  (TM  is  Cu  or  Zn)  have  been 
analyzed  using  neutron  diffraction  and  ab  initio  molecular  dynamics  (QMD)  modeling.  All  six 
partial  radial  distribution  functions  (PRDFs),  gij(r),  have  been  identified  for  each  alloy. 

2.  The  nearest-neighbor  mode  bond  lengths  are  shorter  than  those  in  competing  crystals.  A 
noticeable  shortening  of  Ca-Cu,  Mg-Cu,  and  Ca-Zn  bond  distances  indicate  strong  interactions 
between  these  atom  pairs,  as  supported  by  the  calculated  electronic  structure.  It  is  suggested 
that  the  bond  shortening  is  enabled  by  the  absence  of  long-range  atomic  order,  which  lowers 
the  free  energy  of  metallic  glasses  and  increases  GFA. 

4.  Pronounced  chemical  short  range  ordering  (CSRO)  near  TM  atoms,  chemical  short  range 
clustering  (CSRC)  near  Ca  atoms  and  a  neutral  environment  near  Mg  atoms  are  shown. 
Increasing  the  Cu  concentration  from  15  to  35  atomic  percent  increases  the  total  coordination 
number  around  Ca,  CNca,  from  13.6  to  15.0,  while  CNMg  and  CNcu  remain  unchanged  at  -12.5 
and  -10.4,  respectively.  The  partial  coordination  numbers  depend  nearly  linearly  on  alloy 
composition,  so  that  the  number  of  Ca  atoms  decrease,  the  number  of  Cu  atoms  increase  and 
the  number  of  Mg  atoms  are  almost  constant  with  increasing  Cu  concentration. 

5.  Voronoi  tessellation  show  that  many  types  of  coordination  polyhedra  are  present,  but  the  most 
common  are  Kasper  and  distorted  Kasper  polyhedra  at  the  favorable  CNs.  The  fractions  and 
distributions  of  these  clusters  depend  on  alloy  composition.  Polytetrahedral-type  clusters  and 
five-coordinated  vertices  dominate  in  the  amorphous  structures,  which  indicate  that  tetrahedra 
and  pentagonal  bi-pyramids  are  the  main  building  blocks  in  these  amorphous  alloys. 
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TABLES 


Table  1.  Density  (in  g/cnr  and  atoms/ A3)  of  Ca-Mg-TM  (TM  is  Cu  or  Zn)  amorphous  alloys 
studied  in  this  work. 


Alloy 

Density 

g/cnr 

Atoms/A3 

Ca4oMg25Cu35 

2.936  ±0.007 

0.03987 

Ca5oMg25Cu25 

2.439  ±0.003 

0.03497 

Ca6oMg25Cui5 

2.039  ±0.003 

0.03097 

Ca6oMgi5Cu25 

2.367  ±0.004 

0.03271 

Ca6oMg25Zni5 

2.004  ±0.004 

0.03023 

Ca6oMgi5Zn25 

2.289  ±0.007 

0.03130 

Table  2.  The  most  frequent  (mode)  distances  (in  Angstroms)  between  pair  atoms  in  Ca-Mg-TM 


(TM  =  Cu  or  Zn)  amorphous  alloys,  in  accord  to  QMD  simulations. 


Alloy 

Ca-Ca 

Ca-Mg 

Ca-TM 

Mg-Mg 

Mg-TM 

TM-TM 

Ca60Mgi5Zn25 

3.81 

3.44 

3.14 

3.16 

2.83 

2.67 

Ca6oMg25Zni5 

3.82 

3.48 

3.14 

3.18 

2.94 

2.57 

Ca6oMgi5Cu25 

3.78 

3.43 

3 

3.14 

2.72 

2.54 

Ca60Mg25Cui5 

3.78 

3.46 

3.03 

3.08 

2.72 

2.51 

Ca5oMg25Cu25 

3.75 

3.41 

3.07 

3.13 

2.75 

2.48 

Ca4oMg25Cu35 

3.75 

3.42 

3.02 

3.14 

2.71 

2.47 
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Table  3.  Weighted  average  (mean)  distances  (in  Angstroms)  between  pair  atoms  in  the  first  shell 


in  Ca-Mg-TM  (TM  =  Cu  or  Zn)  amorphous  alloys,  in  accord  to  QMD  simulations. 


Alloy 

Ca-Ca 

Ca-Mg 

Ca-TM 

Mg-Mg 

Mg-TM 

TM-TM 

CaeoMgisZms 

4.00 

3.67 

3.30 

3.27 

3.06 

2.83 

Ca6oMg25Zni5 

4.01 

3.65 

3.31 

3.37 

3.05 

2.76 

Ca60Mgi5Cu25 

4.00 

3.67 

3.18 

3.25 

2.88 

2.60 

Ca60Mg25Cui5 

4.00 

3.64 

3.19 

3.29 

2.87 

2.60 

Ca5oMg25Cu25 

4.01 

3.62 

3.25 

3.24 

2.89 

2.58 

Ca4oMg25Cu35 

3.95 

3.64 

3.21 

3.31 

2.83 

2.59 

Table  4.  The  first  shell  cut-off  distances  between  pair  atoms  (in  Angstroms)  in  Ca-Mg-TM 


amorphous  alloys  (TM  =  Cu  or  Zn),  in  accord  to  QMD  simulations. 


Alloy 

Ca-Ca 

Ca-Mg 

Ca-TM 

Mg-Mg 

Mg-TM 

TM-TM 

Ca6oMgi5Zn25 

4.91 

4.55 

4.21 

4.33 

4.13 

3.46 

Ca6oMg25Zni5 

4.95 

4.55 

4.2 

4.33 

3.85 

3.46 

Ca6oMgi5Cu25 

4.93 

4.72 

4.19 

4.01 

3.84 

3.33 

Ca60Mg25Cui5 

4.94 

4.58 

4.13 

4.05 

3.89 

3.31 

Ca5oMg25CU25 

4.99 

4.58 

4.4 

3.97 

3.84 

3.42 

Ca4oMg25Cu35 

4.88 

4.63 

4.26 

4.22 

3.58 

3.25 
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Table  5.  Metallic,  rm,  [43]  and  covalent,  rc,  [44]  crystalline  bond  distances  between  Ca,  Mg  and 
Cu(Zn)  atom  pairs.  These  distances  are  compared  to  rMode  and  rMean  from  the  QMD  PRDFs.  rMode 
and  rMean  are  the  average  values  for  the  studied  alloys. 


Ca-Ca 

Ca-Mg 

Ca-Cu 

Mg-Mg 

Mg-Cu 

Cu-Cu 

Ca-Zn 

Mg-Zn 

Zn-Zn 

Tm  (A) 

3.94 

3.57 

3.25 

3.20 

2.88 

2.56 

3.31 

2.94 

2.68 

rc  (A) 

3.52 

3.17 

3.08 

2.82 

2.73 

2.64 

2.98 

2.63 

2.44 

I*Mode  C^-) 

3.78 

3.44 

3.03 

3.14 

2.73 

2.50 

3.14 

2.89 

2.62 

I*Mean  (-A.) 

4.00 

3.65 

3.21 

3.29 

2.87 

2.59 

3.31 

3.06 

2.80 

100  /o(rMean/l*Mode"  1 ) 

5.82 

6.10 

5.94 

4.78 

5.13 

3.60 

5.41 

5.88 

6.87 

1 00%(rM/rMode- 1 ) 

4.23 

3.78 

7.26 

1.91 

5.49 

2.40 

5.41 

1.73 

2.29 

1 00%(rc/rMode- 1 ) 

-6.88 

-7.85 

1.65 

-10.19 

0.00 

5.60 

-5.10 

-9.00 

-6.87 

100%(rM/rMean-l) 

-1.50 

-2.19 

1.25 

-2.74 

0.35 

-1.16 

0.00 

-3.92 

-4.29 

Table  6.  Partial  and  total  coordination  numbers  around  Ca,  Mg  and  TM  (TM  =  Cu  or  Zn)  atoms 


in  the  Ca-Mg-Zn  and  Ca-Mg-Cu  amorphous  alloys,  in  accord  to  QMD  simulations. 


CN 

Ca6oMgi5Zn25 

Ca6oMg25Zni5 

Ca6oMgi5Cu25 

Ca60Mg25Cui5 

Ca50Mg25CU25 

Ca40Mg25CU35 

Ca-Ca 

9.3 

9.1 

9.3 

9.2 

8.6 

7.0 

Ca-Mg 

2.1 

3.4 

2.3 

3.5 

3.8 

4.3 

Ca-TM 

3.0 

1.8 

2.8 

1.6 

3.2 

4.9 

Total 

14.4 

14.3 

14.4 

14.3 

15.6 

16.2 

Mg-Ca 

8.5 

8.3 

9.0 

8.5 

7.7 

6.9 

Mg-Mg 

1.5 

2.4 

1.0 

1.8 

2.1 

2.6 

Mg-TM 

1.9 

0.8 

2.0 

1.1 

2.1 

3.2 

Total 

11.9 

11.5 

12.0 

11.3 

11.8 

12.8 

TM-Ca 

7.3 

7.1 

6.7 

6.4 

6.4 

5.6 

TM-Mg 

1.1 

1.4 

1.2 

1.8 

2.1 

2.3 

TM-TM 

0.8 

0.8 

1.0 

0.6 

1.0 

1.6 

Total 

9.2 

9.3 

8.9 

8.8 

9.5 

9.4 
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Table  7.  The  short-range  order  parameters,  oti(jk),  for  the  first  coordination  shell  of  Ca,  Mg  and 
TM  atoms  in  the  Ca-Mg-TM  amorphous  alloys  (TM  is  Cu  or  Zn). 


Ca 

Mg 

TM 

Ca60Mgi5Zn25 

0.11 

-0.03 

-0.22 

Ca60Mg25Zn15 

0.09 

-0.06 

-0.08 

Ca60Mg15Cu25 

0.11 

-008 

-0.18 

Ca60Mg25Cui5 

0.11 

-0.13 

-0.10 

Ca5oMg25Cu25 

0.10 

-0.11 

-0.19 

Ca40Mg25CU35 

0.05 

-0.05 

-0.29 

Table  8.  Fraction  of  clusters  identified  in  the  amorphous  structures  of  Ca-Mg-TM  alloys.  Only 


clusters  populated  with  the  fractions  of  0.1  or  higher  in  at  least  one  of  the  alloys  are  tabulated. 


Cluster  Type 

Ca6oMgi5Zn25 

Ca6oMg25Zni5 

Ca6oMgi5Cu25 

Ca60Mg25Cui5 

Ca50Mg25CU25 

Ca4oMg25Cu35 

Ca-centered  clusters 

(0,2, 8, 4) 

0.12 

0.15 

0.08 

0.10 

0.05 

0.01 

(0,1,10,4) 

0.04 

0.11 

0.18 

0.11 

0.07 

0.04 

(0,2, 8, 5) 

0.04 

0.10 

0.06 

0.04 

0.10 

0.03 

Mg-centered  clusters 

e imac 

0.23 

0.13 

0.00 

mmm 

0.10 

0.13 

0.10 

0.00 

(0,0,12,0) 

0.20 

0.08 

(0,2, 8, 2) 

0.10 

0.18 

0.02 

(0,3, 6, 4) 

0.07 

0.16 

0.17 

0.00 

0.16 

0.06 

(0,1,10,2) 

0.10 

0.06 

0.03 

0.08 

0.04 

0.20 

TM-centered  clusters 

(0,3, 6,0) 

0.30 

0.13 

0.24 

0.27 

0.16 

0.11 

(1,3, 3, 3) 

0.08 

0.00 

0.10 

0.02 

0.06 

0.07 

(0,3,6, 1) 

0.20 

0.10 

0.14 

0.10 

0.10 

0.04 

(0,2, 8,0) 

0.16 

0.27 

0.12 

0.03 

0.12 

0.11 

(0,2,8, 1) 

0.08 

0.20 

0.06 

0.03 

0.04 

0.10 
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Table  9.  Global  and  atomic  packing  fractions  in  Ca-Mg-TM  amorphous  alloys. 


Alloy 

Packing 

fraction 

Average  Voronoi  volume  (A3) 

Atomic  packing  fraction 

Ca 

Mg 

TM 

Ca 

Mg 

TM 

Ca6oMgi5Zn25 

0.76 

41.0 

24.2 

14.3 

0.78 

0.71 

0.71 

Ca6oMg25Zni5 

0.76 

40.6 

25.3 

15.9 

0.79 

0.68 

0.64 

CaeoMgi5Cu25 

0.78 

39.9 

23.7 

12.4 

0.80 

0.72 

0.71 

CaeoMg25Cui5 

0.77 

40.5 

24.2 

12.8 

0.79 

0.71 

0.68 

Ca5oMg25CU25 

0.79 

38.8 

24.1 

12.6 

0.83 

0.71 

0.70 

Ca40Mg25CU35 

0.80 

38.0 

22.6 

12.0 

0.84 

0.76 

0.73 

Table  10.  Cluster  packing  fractions  in  Ca-Mg-TM  amorphous  alloys. 


Alloy 

Average  vertex 
coordination 

Average  cluster  volume 

(A3) 

Cluster  packing  fraction 

Ca- 

Mg- 

TM- 

Ca- 

Mg- 

TM- 

Ca- 

Mg- 

TM- 

Ca6oMgi5Zn25 

5.13 

4.97 

4.78 

107.6 

83.3 

64.0 

0.77 

0.75 

0.76 

Ca6oMg25Zni5 

5.14 

5.00 

4.84 

108.5 

84.0 

65.1 

0.77 

0.74 

0.74 

Ca6oMgi5Cu25 

5.14 

5.00 

4.77 

104.2 

81.9 

57.9 

0.79 

0.77 

0.77 

CaeoMg25Cui5 

5.12 

4.96 

4.77 

107.7 

82.0 

58.7 

0.78 

0.76 

0.76 

Ca50Mg25Cu25 

5.17 

5.02 

4.80 

103.7 

77.6 

58.4 

0.80 

0.77 

0.78 

Ca4oMg25Cu35 

5.16 

5.04 

4.86 

96.9 

73.6 

53.4 

0.82 

0.80 

0.80 
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FIGURE  CAPTIONS 


Figure  1.  Total  neutron  scattering  factors  of  several  Ca-Mg-Zn  and  Ca-Mg-Cu  amorphous  alloys. 
Figure  2.  Experimental  (solid  lines)  and  QMD-simulated  (dashed  lines)  total  radial  distribution 
functions,  G(r),  of  Ca-Mg-Zn  and  Ca-Mg-Cu  amorphous  alloys. 

Figure  3.  QMD-simulated  partial  radial  distribution  functions,  gij(r),  for  (a)  Ca6oMgi5Zn25  and 
(b)  Ca6oMgi5Cu25  amorphous  alloys. 

Figure  4.  Composition  dependence  of  two  PRDFs,  (a)  gcu-Cu(r)  and  (b)  gca-Mg(r)  for  Ca-Mg-Cu 
amorphous  alloys. 

Figure  5.  Site-  and  orbital-projected  electronic  density  of  states  for  (a)  Ca6oMgi5Zn25  and  (b) 
Ca60Mgi5Cu25  amorphous  alloys. 
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1.6 


Ca60Mg25Zni5 


Figure  1.  Total  neutron  scattering  factors  of  several  Ca-Mg-Zn  and  Ca-Mg-Cu  amorphous  alloys. 


Figure  2.  Experimental  (solid  lines)  and  QMD-simulated  (dashed  lines)  total  radial  distribution 
functions,  G(r),  of  Ca-Mg-Zn  and  Ca-Mg-Cu  amorphous  alloys. 
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Figure  3.  QMD-simulated  partial  radial  distribution  functions,  gij(r),  for  (a)  Ca6oMgi5Zn25  and 
(b)  Ca6oMgi5Cu25  amorphous  alloys. 
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Figure  4.  Composition  dependence  of  two  PRDFs,  (a)  gcu-Cu(r)  and  (b)  gca-Mg(r)  for  Ca-Mg-Cu 
amorphous  alloys. 
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Figure  5.  Site-  and  orbital-projected  electronic  density  of  states  for  (a)  Ca6oMgi5Zn25  and  (b) 
Ca6oMgi5Cu25  amorphous  alloys. 
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